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SUMMARY 

Revised  and  newly  dev  eloped  formulas  tor  prediction  ot  plastic  del  orinations  ot  struc- 
tural  elements  in  suppressive  structures  are  presented  in  this  report.  These  formulas  supple- 
ment previously  reported  design  equations.  Revised  equations  give  improved  agreement  with 
literature  data,  while  equations  for  additional  structural  elements  allow  prediction  for  a 
wider  spectrum  ot  structural  components.  Some  comparisons  are  made  ot  predictions  from 
these  design  formulas  with  results  of  computer  code  calculations. 
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ADDITIONAL  ENERGY  SOLUTIONS  FOR  PREDICTING 
STRUCTURAL  DEFORMATIONS 


I INTRODUCTION 

A number  ot  structural  uv-Mg-  . I'..."  ' '»•**••'  •.'••‘••'opeil  >11  the  siijxnr-*\ov  e-t**  ; - 

lures  program  lo  picdict  the  plashe  detormations  ot  \anous  struetural  elements  to  blast 
loads  These  ecjualious  are  based  on  the  principles  ol  limit  design  which  allow  errge  plastic 
deformation  sliorl  ol  collapse  or  failure,  and  on  cncigx  balance  principles.  Reierenccs  I and 
2 report  our  initial  ctloils  m ilexeh'pmeiil  x>l  Mich  siriiclui.il  response  expiations.  and  include 
design  formulas  lor  predicting  permanent  delormations  ol  simple  rigid  piastre  systems, 
beams  vs  it li  xatixms  boundarx  conditions,  circular  plates,  reciangiil.it  plates,  spherical  shells, 
and  cylindrical  shells  In  t lie  earlier  xvork.  we  reported  onlx  expiations  tor  either  asymptote 
of  blast  loading,  i.  c . the  nnpulsixe  loading  regime  or  the  q.uM-static  loading  regime 

This  report  extends  and  supplements  the  work  ot  Reierenccs  I and  2.  It  includes  revised 
analyses  ol  plate  iespoi.se  which  acinexc  better  agreement  with  experimental  data  from  the 
literature,  and  new  analyses  lor  other  structural  elements  winch  luxe  been  eneoi'ni.Te  l in 
the  design  of  suppressne  structures  Some  x'l  the  tormuljs  applx  for  the  simultaneous 
application  of  the  two  asymptotic  loading  conditions 

This  xxork  was  performed  lor  I dgcwooj  Arsenal  under t ontract  l)AAA7>-l ' < as 

part  of  the  suppressive  structures  program 


II.  ANALYSIS 


A.  General 

Most  ol  the  analyses  which  we  luxe  conducted  in  the  suppressive  structures  program  are 
based  on  energy  solutions  to  estimate  the  ultimate  load-carrying  ability  ol  structural  com- 
ponents. The  technique  used  to  create  these  analyses  was  developed  in  References  1 and  2. 

A critical  first  step  in  this  approach  is  to  assume  a final  deformed  shape  of  a structural  element 
so  the  strain  energy  can  be  computed.  In  the  impulsive  loading  realm,  one  then  computes 
the  kinetic  eneigy  imparted  to  the  structure  aiul  equates  tins  to  the  strain  energy.  In  the 
quasi-static  loading  realm,  one  computes  the  wuik  pcifoimeu  by  tin-  peak  force  deflecting 
the  structure  and  equates  this  to  the  strain  energy.  References  1 and  2 used  elementary 
rheological  models,  and  data  for  beams  and  plates  to  illustrate  the  procedures.  This  report  is 
a supplement  cn  the  .espouse  ol  plates  and  other  structural  element-.  Earlier  analyses  for 
beams  are  still  valid;  however,  plate  solutions  should  be  replaced  with  the  following  dis- 
C'-'oii.  as  additional  data  luxe  led  to  modifications  that  give  mote  accurate  numerical 
result. 
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B,  Clumped  Circular  Plates 


Tlie  Inst  plate  solution  will  he  lor  . lamped  eireulai  plates  suhjei  ted  to  either  umlorm 
impulse  or  umlorm  pressiiies  Tins  solution  adds  additional  terms  to  the  strain  energy 
c\p:essions  in  Releren  -i-s  | atul  2.  as  both  bending  and  i xtensional  action  are  present  An 
aj  , 'opnute  assumed  deUirnted  shj|-e  loi  a .lamped  deformed  eiuular  j late  is 

u=t'«'  (‘  a>)  iU 


Idle  defonned  shape  is  desenbed  b>  a polar  eoordmate  system  with  its  origin  at  the  center  ol 
the  plate  Bsvausc  ol  symnv.rx.  the  delonned  slupe  is  independent  ot  the  angle  0 and  no 
shears  exist  Iquation  ( I ) ;s  an  equation  lor  a eone  it  is  an  acceptable  detormed  shape,  as 
at  r - 0.  the  maximum  deformation  is  finite  tie  - »»,,  >.  ant  at  r = R . » = 0.  The  slope  is  not 
zero  at  r - R.  as  a plastically  ik  forming  plate  is  bent  ovci  tile  edge  ol  the  clamp  along  an 
abrupt  y ield  line  that  lias  an  angle  change  ol  4»  = U (d>.  dm,  = A . or  - (wt,  j*\  » 


Because  no  change  in  length  occurs  circttmlerenti  ill;’  at  the  boundary . there  is  no  cir- 
cumferential extensional  strain . Itowc.cr,  circumferential  bending  strain  docs  occur  The 
circumferential  plastic  strain  energy  per  unit  differential  length  dr  equals  the  plastic  yield 
moment  limes  the  circumference  tunes  the  circumferential  curs jturc  or  to,  lr  4)  times 
( 2 trn  times  | ( 1 ritdu  dr)j . Radially,  the  deforming  plate  stores  energy  in  membrane  action, 
litis  lauiai  suam  cneigs  pet  unit  uilicicntiai  length  dr  equals  iue  yield  s'uess  times  the 
thickness  times  the  circumference  of  a ting  at  location  r.  times  the  extensional  strain.  Using 
the  well-known  approximate  extensional  strain  relationship  1 2 tdu  dr)J  lor  the  average 
change  in  radial  strain  means  that  the  radial  extensional  strain  energy  stored  in  the  plate  per 
unit  differential  length  dr  is  given  by  to,  i tunes  1 2rrr/t ) times  1 1 /2 1 (dn\  dr): . The  third  and 
final  amount  ol  strain  energy  stored  in  a clamped  plate  is  that  associated  with  plastically 
bending  the  plate  user  the  edge  of  the  clamp.  Bending  strjin  energy  along  this  yield  line  is 
vdlcuiated  by  i.ivdiipiy  mg  the  y ield  moment  per  unit  width  times  the  circumference  of  the 
plate  times  the  change  in  angle  at  the  edge  of  the  plate  ot  (o,  lr  ,4 ) times  ( 2ff/\ ) times  ln„  R ) 
The  total  strain  energy  i is  the  sum  of  these  contributions  or 


U=  )(>«,  (”*')  4 J-(2rrr,  ^ +J  »«,  »<2*r/,)  (’)(*-)  dr 


bending  over 
tint 


circumferential  bending 


radtal  extension 


(2) 


I'.ffcrcntiatmg  Lq.  ( I ).  substituting  it  into  f.q  ( 2 ).  and  gathering  terms  men  yields. 


, . * , , rr  * dr  , , * rdr 

I = 0-  lr  n..  + o h- u..  f + jru,/tn'  t 

2 ' 2 - R - R2 

o O 


(1) 


8 


Or: 


H 

('  - K4I,./|:  n„  -t  , 0,  /m,', 


<4) 


The  k t ne  tic  energy  A h to:  a uniformly  applied  impulsive  load  imparted  to  a plate  in  ob- 
tained by  summing  the  impulse  squared  divided  b>  two  times  the  incremental  mans  over  the 
surface  ot  the  entire  plate.  This  procedure  leads  to  the  tollowtng  integration. 


* i } (’trr)J  (dr): 

it  = r _ _ _ 

i VK^hdr) 


Or 


kl.  = 


*/;  a,: 

7p/T 


<5) 


(0) 


b qua  ling  <•  ( bq  (4  >]  to  Af  | bq.  (<>))  y idds  the  asv  inptote  lor  the  impulsive  loading 
realm. 


’t"|  ; ) + ( / ) | circular  plate,  impulsive  realm  | (7) 

A comparison  between  b.q  < 7 > (ora  uniform  impulse  imparted  to  a clamped  circular  plate 
and  experimental  t est  Jala  can  be  made  using  results  by  F lorence.’  Kesidual  pennanent  mid- 
span deformations  were  measured  on  clamped  circular  6061-T6  aluminum  plates  and  1018- 
cold  rolled  steel  plates  that  had  been  load  'd  uniformly  with  various  thicknesses  of  sheet 
explosive.  1 he  7 2 aluminum  data  points.  20  steel  data  points,  and  bq  { 7 ) are  all  show  n in 
Figure  I There  is  a reasonable  agreement  between  bq  and  the  test  results.  It  there 
ex. .its  a systematic  error,  it  is  a tendency  for  the  analytical  curve  to  slightly  underestimate 
defc. nauons  whenever  »*„  7/  is  large  This  error  is  probably  caused  by  the  assumed  deformed 
shape  not  yielding  the  minimum  strain  energy  bor  small  values  o.  ><  „ ti.  the  analy  ticjl  curve 
overestimates detormations  Thu  error  is  probably  created  because  we  assume  that  defor- 
mations extend  over  the  entire  span  ot  the  plate  When  loads  and  deformations  are  small 
the  deformed  shape  undoubtedly  covers  only  a portion  of  the  entire  plate. 

A solution  for  a clamped  circular  plate  in  the  quasi-statu  lo.  dmg  realm  is  obtained  by 
computing  the  work  and  equating  it  to  the  strain  energy . The  work  It  7.  is  given  by  the 
integral  of  the  maximum  pressure  times  a differential  area  times  the  deflection  of  the 
differential  area,  and  is  expressed  by  bq  (,'if 


n 

hx  - /V  ’irrdr 

O 


(8) 
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pn  wilt  a anil  i.-  v.'ati.laill.ii  Sim  il  a t j u:i  by  ill  ail  a ti'Hip  noiii.  I a -t  K ->nl  !•>  b>  / illia>  11  I’lii  illai>  1 . 
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/-  DELTA  FUNCT I OfJ  ( Ignored  as  insignificant  energy 
2280  psi  -j  / input  relative  to  that  from  a quasi -static  loading  > 

I ( msec  ) _| 

PEAK  OIJASI-STATIC  PRESSURE 

•-CALCULATED  AVERAGE 
' QUASI-STATIC  PRESSURE 


I li .1  HI  ; lit  AM  Kl  11  I UMUSi.  IN  V.Kia  i IKt  l l \K  Pi  V 11  I I SI  -> 

l‘io|  I ed  I'l  1 Ijtnle  ,tl  V1  il.lt.  1 poult'-  I tom  Kek'.eik  4 t;lit  t t{  i I 2 I A)Ti.'i'i:M1I  l\  c\ 
eeilent  o\ei  1 1 u'  entile  r.mjte  Snhse»iueni  ted  results  on  senate  an>i  n\  tjnguiai  piates.  iuadeu 
piiemnaiiealh  so  no  i e\ «.  ’her  a lions  or  n Heeled  sIuk  ks  ate  i omul,  will  show  that  no  I ' 
tael-  i !'  rotjmrotl  to;  overshoot. 

C Ciamped  Rectangular  Plates 

1 lie  hum  general  plate  solution  is  loi  the  pladu  response  ol  lectangular  ( l.<’e>  t j either 
unilormU  applied  inipnl.es  or  pressures  This  solution  i:r.tod:ues  new  s omplk  aliotts  t's.ause 
ol  esiension.il  ami  hemline  shea:',  l ot  a v.  Ltnped  u\  tanituiai  plnU  . hem!  mg  w II  o»  em  alone 
>ieU  I iiu  e at  the  I oniulattes.  and  the  plate  will  delo.m  in  memhtaiu  action.  At.  apptopuate 
ile tor  lied  shape  is  assumed  to  he  given  h\ 


n 


it  \ 


:\ 


ri 


< i ,'t 


u here 


A and  1 at  e hall  span* 

v aiul  i ate  j rectangular  cooiJiiuie  > \ > t e 1 1 * with  origin  ai  'lie 
eenier  ol  the  plate 

1 Ins  assiniu  J Je  lor  med  dupe  meets  the  appiopr  i.ue  d elles  lion  „ i iter  i.-  in  t hi  in  id'll  -e  and 
along  the  homul.iues  J he  slope  in  the  muld!.-  is  /er-'.  -./Ink-  ihe  dope  at  the  ho.Jiklat.es  is 
ilk  negative  oi  the  jni;le  ehanjre  when  the  plate  is  hen:  met  the  ei.nnp.  T he  sti  Jill  cnergv 


ir 


per  unit  volume  in  a structural  element  under  a biaxial  state  ol  stress  is 
l 


Vol. 


/ [°x  > dtj,  ( 4-  2oxl  dt  rl  4 u , , dc,  ,] 


(14) 


Because  we  have  yielding,  we  will  assume  oxx  - o,  and  c,  = u,  . bu:  lor  the  shearing  stress 
we  will  use  a Huber-Mises-Hencky  distortion  energy  yield  criteria  ol  o(l  = (a,  />/T).  At  lirst 
glance,  one  would  say  that  Mohr’s  circle  for  stress  is  a point  it  o, , = o,  x ; hence,  no  shears 
can  occur.  Such  an  observation  is  correct  lor  any  instant  in  time.  but  normal  stress  dis- 
turbances propagate  faster  than  shearing  disturbances.  We  are  writing  an  energy  expression 
for  final  state  and  do  not  care  about  the  timing  of  events  it  is  sufficient  for  our  purposes 
to  know  that  plate  distortions  from  normal  stresses  and  lrom  shearing  stiesses  can  occur  at 
different  times;  thus,  permitting  the  use  of  distortion  energy  principles  to  determine  yield 
stresses.  The  normal  extensional  strains  in  a rectangular  plate  arc  similar  to  those  in  a circular 
plate  with  cvv  = ( 1 ’)liw  i)u;  and  evv  = < 12)(dn  , di  )J . The  extensional  shearing  strain 
t v , = (Die/ d.v)(du  /t).\ ) as  a first  approximation  These  observations  mean  thal  the  exten- 
sions I strain  energy  L 'e  is  given  hv : 


X 1 h 

- 4 ('  d v f d i |‘  dc  « 

t>  o o 


' 1 - / \ /dw\/du\ 

+ 4 f dv  f dr  I dr  < 2 ){  -L)  ( - (15) 

\s/4/  \&X  l\iU  / 


u o t> 


DM lerentiating  the  deformed  shape  | bq.  ( 14)1 , substituting  it  into  hq  (15),  and  performing 
the  required  triple  integration  y ields: 


<■  ~ «,  /'>>;, 
O 


Y 


Al  4 

4-  o, 


/in;, 


(lb) 


Normal  stresses  associated  with  bending  of  the  yield  line  along  the  clamp  cause  additional 
strain  energy  to  he  stored  1 his  energy,  analogous  to  similar  terms  lor  the  circular  plate, 
is  given  by : 


dv 


<I7) 


Substi'uting  <o,  ) 4 for .'/,  . differentiating  I he  deformed  shape  so  it  can  be  evaluated  at 

the  boundary  and  further  substituting  then  yields 
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(18) 


ov  /i3 


JTV 

cos  dv 


Or: 


i/ 


hll 


<ll/> 


Finally,  there  exists  a shearing  strain  energy  associated  witli  bending  £//,,  in  the  localized 
zone  around  the  yield  lines.  II  the  small  deformation  equation  lor  shear  strain,  e, , ■ - 
2 : (0:  u 'dvd.i ),  is  used,  the  energy  i f, s per  unit  volume  is  given  by : 


I'h, 

Vol 


■to  (■" 


d : »e  \ 
3 v 3 1 / 


(20) 


at  the  boundaries 


If  we  assume  that  the  localized  shear  zane  is  four  plate  thicknesses  wide,  Eq.  (20)  becomes 

4 hi  2 > 


L'h  = 


/ d:  n \ 

-S-Or  (■».»  | 2d:  f d.l  (4/1)0  --e 

a J \Zx*y/x . v 

4 o,  (4)  f 2d:  f d v (4/1)0  . - - ) 


(21) 


Or,  alter  differentiating  Eq.  (13)  to  substitute  for  d:u73v3r  at  the  boundary 


f-  hi  = 


3 27T-  0,-llW, 


J7  XY 


hl2  l r mri  " - ' r irvri 

/ J - sin  — d:  d.r  x f f : sin  — | d’  dv 

i .1  L —'-sJ 


h'2  A 


(22) 


Performing  the  required  double  integration  then  yields 


I tor  , I"  /,•  /i~| 

L'ht  \/5  0,/'  Uu  b + > J 


(23) 


Idle  total  strain  energy  ( stored  in  a rectangular  plate  is  the  sum  of  i'r  plus  plus  r<,s  or: 

(24) 
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The  kinetic  energy  Kh  imparted  to  a plate  is  not  dependent  upon  the  deformed  shape. 
Eor  a uniformly  applied  impulsive  load  imparted  to  a rectangular  plate,  the  kinetic  energy 


15 


is  obtained  by  summing  up  the  impulse  squared  divided  by  two  times  the  incremental  mass 
over  the  surface  of  the  plate.  The  appropriate  integration  is  given  by: 


x > 


**  = 4 f / 

<>  O 


i;  (d.v)2  (di ): 
2 ph  (d.v)  (d.v) 


or 


A E = 


2 <;  X Y 

ph 


(25) 


<26) 


tquatmg  [ tq.  (24))  to  A A [tq.  (26))  y ields  the  asymptote  for  the  impulsive  loading 
realm. 


fiiii  l =i(^L(>  yl  + ^( -)[-]  [i  j i'f 

pof/i«  2 V //  / j_  V.v/J  s/T  \ h / |_vj  L *J  16  V/!  / _ V .v/J 


4 


(27) 


Tins  is  the  geueial  solution  for  an  impulsively  loaded  plate.  It  can  be  compared  to  test 
data  reported  by  Jones,  et  al.s  Rectangular  plates  with  an  aspect  ratio  Y/X  equal  to  1.695 
were  loaded  with  sheet  explosive  in  these  tests.  Both  hot-rolled  mild  steel  plates  and  6061-T6 
aluminum  plates  were  iested  and  can  be  seen  plotted  in  Figure  4.  Although  various  thick- 
nesses of  plate  were  used,  we  will  use  an  average  value  for  (Xlh)dy),  = 10.2  in  the  second  term 
of  tq.  (27).  Using  an  average  value  for  X/h  creates  very  little  error,  as  the  (w0!h)i  membrane 
term  predominates.  Substituting  for.'.’/A  in  tq.  ( 27 ) y ields  tq  (28),  which  is  compared  to 
the  test  data  in  Figure  4 


5.76 


(28) 


Agreement  is  excellent  when  data  are  compared  to  tq.  (28)  as  in  Figure  4.  The  little  error 
which  does  oc  ur  for  small  values  of  »»■„//;  is  piobably  caused  by  the  deformed  shape 
occupying  only  a portion  of  the  full  span  when  deformations  are  small. 


A quasi-static  loading  realm  solution  for  rectangular  plates  is  developed  if  the  work  is 
equated  to  strain  energy  The  work  WA  is  given  by  tq.  I 29)  in  a manner  analogous  to  that 
for  a circular  plate. 


v > 

Wk  = f f 4 Ah  dv  d.v 

O O 


(29) 
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Substituting  the  assumed  deformed  shape  as  given  by  Eq.  (13)  and  performing  the  required 
double  integration  then  yields: 


16 

H'A  =-.-/*  u\,  A'  >' 


(30) 


Equating  U (Eq.  (24))  to  Wk  (Eq.  (30)1  yields  the  asymptote  for  the  quasi-static 
loading  realm. 


4v/T 


We  can  evaluate  the  validity  of  Eq.  (31 ) by  comparing  it  to  test  data  taken  by  Hooke 
and  Rawlings.6  In  these  experiments,  damped  rectangular,  mild-steel  plates  of  various 
thicknesses  and  aspect  (Y/X)  rr.tios  were  subjected  to  step  pneumatic  pressures  of  very  long 
duration.  Because  the  pressures  are  pneumatically  applied,  no  overshoot  is  present  as  in  the 
quasi-statically,  shock-loaded,  circular  plates.  Figures  5 and  6 present  some  of  these  data  in 
plots  of  scaled  applied  load  \PX2  l(ayh2)  1 as  a function  of  scaled  permanent  mid-span 
deflection  (u0//i),  for  square  plates  and  for  rectangular  plates  with  an  ( X\Y ) of  1/2.  An 
average  value  for  scaled  plate  thickness  (A '//;)  of  53.5  was  used  to  determine  the  analytical 
lines  in  Figures  5 and  6.  Because  {X/h)  varied  over  a limited  range,  and  because  its  influence 
is  only  moderately  important,  a single  value  suffices  when  tq.  (31 1 is  applied.  Substituting 
for  (X/Y)  and  (X/h)  in  Eq.  (3 1 ) yields  Eq.  (32)  for  the  square  plate 


PX 2 

Oy  h\ 

and  Eq.  (33)  for  tests  on  a rectangular  plate  with  A'/Fequal  to  1/2. 

= 1.038  + 1.664  ^ y'j  rectangular  plate  tests  by  Hooke  (33) 

Agreement  is  excellent  in  Figures  5 and  6.  Because  test  results  are  now  consistent  with 
the  theoretical  approach  in  both  the  impulsive  anu  quasi-static  loading  realms,  the  physical 
viewpoints  being  reflected  in  this  analysis  appe^*-  m be  substantiated,  and  we  now  feel  we 
know  how  to  handle  structural  members  which  are  under  a biaxial  slate  of  stress  as  well  as  a 
uniaxial  one. 


= 1 .95  + 2.96 


(t) 


square  plate  test  by  Hooke 


(32) 
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D.  Rings  Constraining  I beams 


l,i  t ho  final  design  concept  lor  the  Category  I structure,  the  basic  structure  consists  ol  a 
series  ot  interlocking  vertical  I-beams,  supported  by  several  rings  or  hoops.  Because  ol  this 
symmetry  of  the  structure  and  the  applied  loading,  ring  rcspoc.se  is  represented  by  a uniform 
radial  expansion  It  the  energy  absorbed  in  deformation  of  the  I-beams  is  neglected,  then  the 
ring  response  can  be  conveniently  determined  by  an  energy  balance.  1 o do  this  the  kinetic 
c lergy  imparted  to  the  I-beams  and  nngs  by  the  initial  blast  wave,  plus  the  work  done  by  the 
cp.  isi-stalic  pressure  as  the  cy  linder  expands,  are  equated  to  the  strain  energy  in  the  rings. 

Figure  7 shows  the  geometry  ol  the  cylinder  wall. 

For  uniform  radial  expansion  the  maximum  possible  kinetic  energy  imparted  to  the  beams 
i"vj  ring  is 


.r  "('drV/t*,.. 

vlf  - o 2 \ III  / 


id  _ IK*!?}/ 

Ill 


(35) 


where  in  = mass  per  unit  circumterenee 


- Arp 


Ls'»b 

CSBr, 


Ln  length  ot  the  beams  supported  by  one  ring 

mu  mass  per  unit  length  of  the  beams 

CSBr  circumferential  spacing  ot  the  beams  as  measured  at  R& 
Rr  radius  to  the  ring  (£ 

Ar  - ring  cross-sectional  area 
p - density  of  the  ring  material 
/,  - reflected  specific  impulse 


The  maximum  possible  work  done  by  the  quasi-static  pressure  during  the  ring  expansion 
is  given  by: 


Wk 


?7T 

$ PLB  - Mi  . d C=  f 

<■  O 


R d 0 


2*Rk 


<3(>) 


where 


R - mean  radius  of  the  mteilocking  I-beams 
P - peak  quasi-static  pressure 
A R - radial  expansion  of  the  ring 
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Tins  expression  assumes  a constant  qu.isi-statu  pressure  -:.u  mm  ,..  'onservative  t s> r rings  with 
lout:  response  tunes  The  degree  ol  conservatism  can  he  estimated  (torn  the  response  cals  il- 
lations reported  in  Reference  7.  The  calculated  response  time  tor  a img  in  the  Category  I 
prototype-  (with  a cross-sectional  area  of  110  iir  I is  found  Irom  Reference  7 to  he  about  40 
ills.  Figure  N shows  the  relationship  between  the  response  time  ol  the  ring  and  the  deux  ol 
the  quasi-static  pressure.  An  estimate  ol  the  error  is  made  b>  comparing  the  actual  area 
under  the  pressure-time  curve  (impulse  produced  h\  the  quasi-statu  pressure!  to  the  atea  lor 
constant  /*.  The  dilterence  m areas,  shown  shaded,  approximates  the  error  It  is 

PiAQ)  \(P  + Pt„  t (40) 

T.rror  = + X 100  = 4 'A?r>  , 


HU  RLS  Sr  111  M X I K Ol  I KKOK  INCl’KKt.l)  HX  ASSl'MPIION 
Ol  «>NS1  \M  ot  XSI  -St  AIK  PRISSt  KI  . 


For  rigid-plastic  behavior,  the  strain  r,  will  equal  AA  A’ a and  the  strain  energx  m th 
nng  undergoing  radial  expansion  is  given  by: 


f = .4 A-  f O,  c.  A’ a AO  = > Aa  ,1a 


where  o,  - yield  stre-s  of  the  ring  material 
Now.  by  equating 


A A -f  It  A = 1 


the  following  relationship  is  obtained 


'Ilk 

m{Af(Ot  —Rl.gP)  vA’#,  f 


(38) 


t o estimate  the  degree  of  conservation  in  Eq.  (38)  caused  by  neglecting  the  energy 
absorbed  by  plastic  deformation  of  the  I-beams,  we  again  used  the  response  calculations  for 
the  beams  and  nngs  described  in  Reference  7 Based  on  the  calculated  deformations  for  the 
beams  and  nng  in  the  Category  I prototype  shield,  the  energy  absorbed  by  each  was  com- 
puted tor  comparison  For  the  b’8  X b7  beams  supported  by  nngs  of  cross-section  140  in:, 
the  plastic  deflection  of  the  inner  beams  between  the  center  and  lower  nngs  was  found  to 
be  about  I .6  inches  which  gives 

H',,7  = I o/l  20  = 0.01 33 


An  estimate  of  the  strain  energy  absorbed  by  the  beams  with  this  venter  deflection  is  found 
from  equation  VC.2I ) of  Reference  I 


1 6 My  W0 
L 


\ 


(39) 


where 

Mt  - plastic  moment  of  the  beam  86  2.613  X 1 0fc  m-lb 
.V  - number  of  beams  in  the  inner  layer  = 157 

Therefore  the  energy  absorbed  by  the  inner  beams  was 

t;=  16(2.013  X 10‘  in-lb)  t0  0133)  1157)  = 87.3  X 10*  m-lb 

If  we  conservatively  assume  that  the  outer  layer  of  beams  absorbs  equal  strain  energy  the  total 
energy  absorbed  by  the  beams  is 

i = 2(87.3  X 10-  in-lb)  - 1 74.6  X 10''  in  ib 

The  energy  absorbed  by  the  nngs  is  estimated  from  the  total  radial  deflection  o!  the  center 
ring  from  Reference  7 and  Equation  (37)  From  the  table  A/f  = 18.2  in.  so  that 


bJR 

Rr 


18.2 

283.5 


= 0 0642 
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with 

Ak  I 40  ur 

o,  40,000  psi 

/ AA'\ 

< - 2»  A\-  I*  u ( ) - <i40..'  X I O'-  in-lb 

\ A ' 

I here  loro  1 1 ic  fraction  ot  the  total  energv  absorbeJ  b\  0k-  beams  is 

174  0 X 100 

•;,/  = — --  :i  4 

<••40.-1  i 174.0 

Although  tbo  number  o high.  bceause  (be  assumption  wj'  nude  tli.d  the  outer  la\  et  •>!  beam 
jbsoibs  the  same  amount  ot  otiotgv  .is  the  inner  beams.  it  still  uulu.iles  tlu,  r e is  j go.>J 
deal  ol  eoiiseiN jtivir.  in  l q ( IX)  I von  so.  the  equation  ejii  be  used  tt>  eibtai  .,,iuL  con- 
servative estmutes  ol  the  rule  requirements. 

E.  Hemispheres  (Domes) 

l!  umlonn  loading  is  assumed  on  the  inner  surlaee  ol  the  sphere,  and  it  the  radial  de- 
flection at  the  base  is  assumed  to  be'  equal  to  that  in  the  remainder  ot  the  hemisphere,  then  a 
undorm  radial  expansion  e\aet!>  represents  »he  dome  response  ! Ins  i-  shown  in  f igure 
In  actual  structures  the  hemisphere  usually  will  he  secured  at  its  base  to  a cylindrical  section 
whose  radial  expansion  nu>  not  match  ovactlv  that  of  the  hemisphere.  In  other  eases  the 


dome  ma>  even  be  sccureJ  “rigidly”  to  a iviidorccd  concrete  toundation.  The  fllect  ol  this 
re  ill  anil  on  the  total  energy  absorbed  by  the  dome  should  be  acceptably  a. .all  lo'r  the  “rigid 
c..:,c  and  negligible  foi  the  other.  The  assumption  ol  Uniterm  loul.ng  is  also  acceplab'.  because 
the  quas-static  pressure  is  uniform  „nd  the  initial  blast  wave  is  appioxunatcly  umfortn 
depending  on  the  charge  kx.  lion  relative  to  the  center  of  the  hemisphere  as  well  as  the  charge 
geometry.  We  have  also  observed  that  in  the  C ategory  I design,7  the  qujsi-stutic  prev.uie 
makes  the  greatest  contribution  to  the  dome  thickness  requirements 


As  m the  case  of  rings,  the  maximum  dome  deformation  is  obtained  by  c puling  the 
kinetic  energy  imparted  to  the  dome  by  the  initial  blast  wave,  plus  the  weirk  done  by  the 
ijujsi-stjiic  pressure  during  hemispherical  dome  expansion,  to  tin-  strain  energy  absorbed  by 
the  dome  kinetic  eixigy  is  given  by 


A A = 


t - ml  1 iL-t  _ j — d.1  - 

A - A ‘ 


xR : /; 
m 


(40) 


w hs  re 

>n  mass  per  unit  area  of  tne  dome 
A’  dome  radius 
i,  specific  reflected  impulse 


lhe  work  done  by  the  quasi-static  pressure  dunng  expansion  of  the  dome  is 

If  A = f P AK  d.-l  - lx  A’ * P ( — ) i4 1 I 

.a  V H ' 

R is  the  peak  quasi-static  pressure  and  its  decay  over  the  response  tim  • of  the  dome  has  bernt 
ignored 

l.lastic  strain  energy  per  unit  volume  ol  the  dome  material  isgivui  by  the  well-known 
formula 


I 

■ii  =;loJt1  -re,  c,  + tiy  ) dvd.i  dr 


t4.’j 


bor  rigid,  perlectly-plastic  behavior,  shear  stresses  vanish  and  the  stress  is  constant  so  that 
the  slram  energy  doubles,  fcqujlion  (42)  sail  then  b>'  written 


/ - I (Oj,  e , + o,  t , ) 'i  d.  t 
.a 


(4c  i 


where  h = thickness  of  the  material.  Because  the  stresses  and  strains  are  uniform  throughout 
the  hemisphere  for  the  assumption  we  have  made,  the  equation  further  reduces  to 


U - h J 2 o e d A = 4 ir  R7 h a t 

•t 

Since  the  s',  ain  can  be  expressed  as 


e - 


A R 
R 


(44) 


the  strain  energy  finally  becomes 


V = 4 7T  R 2 


(45) 


for  a dome  thickness  h and  yield  stress  ov . 


Now,  equating. 


KE  + Wk = L 

The  following  relationship  is  obtained 


l'r 

phOhuy  - PR) 


wltere  p/i  ■ as  replaced  m in  the  kinetic  energy  term. 


(4b) 


For  some  applications  to  suppressive  shields,  two  cu  *nc  domes  are  used  with  a 
filler  in  between.  The  purpose  of  this  arrangement  is  to  detea.  : ragmen ts  with  the  inner 
dome  and  filler,  allowing  the  undamaged  outer  dome  to  resist  the  internal  loads.  This  dome 
arrangement  is  given  in  Figure  10.  Roponse  of  this  dome  arrangement  differs  from  that  of 
the  single  dome  in  the  following  ways. 


( 1 l 1 he  response  is  altered  because  of  the  added  mass 

(2)  Hie  inner  dome  stretches  along  with  the  outer  dome 

(3)  The  radial  load  is  reduced  because  of  the  reduction  in  the  loaded  radius, 
blit  the  total  vertical  load  is  unchanged  because  of  the  pressure  acting 
on  the  ring  which  closes  the  bottom  of  the  space  between  the  domes. 

Although  the  inner  dome  expands  along  with  the  outer  dome  and  thus  absorbs  energy  also, 
it  is  also  possible  that  catastrophic  rupture  will  initiate  at  fragment  damage  locations  during 


OUTER  DOME 


VERTICAL 
REACTION  (TYP.) 


FIGURE  10.  GEOMETRY  OF  THE  DOUBLE  DOME  SEPARATED  BY  FILLER  MATERIAL 

its  expansion.  Thus,  the  contribution  of  the  inner  dome  to  the  strain  energy  is  ignored.  Also, 
the  reduction  in  the  loaded  radius  is  ignored  because  the  vertical  loading  is  unchanged.  Thus, 
a conservative  estimate  f the  response  is  obtained  by  modifying  only  the  mass  in  the  kinetic 
energy  term  given  by  £q.  (40).  Instead  of  m = ph  for  a single  dome  we  h .ve 

m = ph  + ma 

where  >nB  ~ added  mass  produced  by  the  inner  dome  and  tiller  material  divided  by  the  area  of 
the  outer  dome. 


Equation  (4(>)  becomes  simply 

if 


(ph  + ma  ) 1 2 h oy  — PR 


--’(f) 


(47) 


for  the  double-donic*  arrangement.  To  date  no  comparisons  with  experiment  have  been  made 
for  either  £tj.  <4(i)  or  (47). 


III.  RESULTS 


The  results  of  the  energy  balance  analyses  presented  here  update  and  extend  those  given 
in  References  I and  2.  Because  these  formulas  now  are  rather  numerous,  we  nave  summarized 


them  in  an  extensive  tabic  (Table  1 ),  giving  a description  of  the  structural  element,  the  assumed 
deformed  shape,  and  the  solution  for  combined  impulsive  and  quasi-static  blast  loads.  These 
formulas  reduce  to  the  quusi-static  and  impulsive  loading  asymptotes  for  i,  = 0 and  p = 0. 
respectively.  Symbols  in  Table  1 are  defined  in  Table  2.  No  dimensions  are  given  because 
the  equations  are  in  essentially  dimensionless  form,  and  correct  for  any  sclf-consistent  set 
of  units 


IV.  DISCUSSION 

The  energy  solutions  presented  in  this  report  should  give  designers  additional  tools  for 
rapid  estimation  of  plastic  deformations  in  suppressive  structures  or  elements  of  these  structures, 
for  the  initial  blast  loading,  for  the  longer  duration  quasi-static  pressures  and  for  both  of  these 
loads  applied  simultaneously.  The  revised  equations  for  plate  deformations  consider  strain 
energy  terms  not  included  in  previous  analyses,  and  also  utilize  deformed  shapes  which  are 
closer  to  experimentally  observed  final  deformations. 

This  report  describes  the  current  status  of  the  work  on  a portion  of  the  suppressive 
structures  program  The  energy  solutions  given  here  will  undoubtedly  be  supplemented  by 
other  solutions  or  modified  as  more  data  on  dynamic  plastic  responses  in  this  program: 
notably  BRL.  NSWC  White  Oak,  and  Corps  of  Engineers  at  Huntsville,  Alabama.  Results  on 
response  prediction  methods  from  these  studies  and  other  work  in  the  literature  should  be 
compared  as  the  program  progresses.  The  energy  methods  are  powerful  and  potentially 
qi.  t?  "wfi.t  be  "iwa  they  yield  relatively  simple  design  formulas.  They  do.  however,  involve 
a number  of  simplifying  assumptions  and  should  not  be  expected  to  yield  the  accuracy  in 
response  prediction  possible  with  a number  of  available  dynamic  response  computer  codes, 
nor  can  they  predict  time  histories  of  deflections  or  strains.  One  should  probably  rely  on 
the  energy  solutions  for  initial  analysis  and  design,  and  to  compare  feasibility  of  different 
concepts.  Then,  designs  can  be  refined  by  specific  eak  illations  using  more  complex  and 
more  nearly  exact  dynamic  response  computer  programs.  This  was  the  procedure  followed 
in  the  design  of  the  Category  1 structure  reported  in  Reference  7. 
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Symbol 

A 

b 

CSRp 

h 

>r 

L 

Lb 

fTIfi 

mB 

MP 

X 

P 

Py 

r 

R 

*R 
A R 

H* 

JC 

X 

V 

Y 
P 

Oy 


TABLE  2.  DEFINITION  OF  SYMBOLS  USED  IN  TABLE  l 

Definition 

beam  cross-sectional  area 
ring  cross-section;*!  area 
loaded  width  of  beam 

circumferential  beam  spacing  in  the  1-beam  cylinder  measured  at  Rr 
thickness  of  plate  or  dome 

specific  reflected  impulse  from  initial  blast  wave,  plus  reflections  if  applicable 

length  of  beam  for  which  the  deformation  is  being  determined 

length  of  beam  which  is  restrained  by  a single  ring  *n  the  1-bean,  cylinder 

mass  per  unit  area  of  the  inner  dome  and  filler  material  for  the  double-dome  roof 

mass  per  unit  length  for  the  beams  in  the  1-beam  cylinder 

beam  plastic  moment 

factor  in  the  beam  equation;. V = 1 for  simple  support,. V = 2 for  damped  support 

quasi-static  pressure 

axial  yield  force  of  the  beam 

radius  to  arbitrary  point  on  a circular  plate 

mean  radius  of  the  hemisphere  (dome) 

mean  radius  of  the  ting  me  the  I-beam  cylinder 

radial  expansion  of  the  ring  or  dome 

lateral  d-  tlect.on  of  a beam  or  plate  at  point  x or  r,  respectively 
center  dcflec.ion  of  a beam  or  plate 

distance  along  the  beam  or  plate,  normally  measured  from  the  center 
short-semi-span  of  iV  plate 

distance  along  plate  center  line  normally  measured  from  the  plate  center 

long  semi-span  of  the  plate 

material  density 

yield  stiength  of  the  material 
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